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Abstract. We present terminal velocities of M33 B-supergiants, obtained from STIS HST spectra as part of our programme to 
investigate the Wind Momentum - Luminosity Relationship (WLR) in the Local Group. Terminal velocities are derived from 
their N V, C IV, and Si IV resonance lines in UV spectra. Comparing with lUE spectra of Galactic B-supergiants we found 
evidence of low metallicity in three of our objects. The terminal velocities are consistent with the corresponding values of 
Galactic stars, except for B-133. For this star we find a very large Voo and a red Si IV component deeper than the blue one, 
that might be an indication of binarity. The average ratio between terminal and turbulent wind velocities is 0.25, well above the 
value found for Galactic stars. 
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1. Introduction 

The Wind Momentum-Luminosity 
Kudritzki et al. 1995|; [Puis et al. 1996| , 



Relationship (WLR, 
for recent reviews see 
Kudritzki 1998 and Kudritzki & Puis 2000) offers the very 
attractive opportunity to derive the stellar luminosity directly 
from the analysis of the observed spectrum, provided that the 
wind is radiatively driven. Accurate individual stellar distances 
can then be obtained if the apparent magnitude and extinction 
are known. This method, that has been extensively described 
elsewhere (see for example the reviews cited above), can be 
potentially applied to individual stars up to distances of 1 0-20 
Mpc, reaching the Virgo and Fornax clusters with accuracies 
as low as 0. 1 mag, using present day telescop es and techniques 
([McCarthy et al. 1997|; 



in dif ferent environments w i thin the Local Group (Puis et al. 
1996; [VlcCarthy et al. 1997|; [Kudritzki et al. 1999|; Herrero et 



Kudritzki et al. 1999) 



However, accurate application of the method requires pre- 
vious calibrating work, as the W LR depends on the stellar 
metallicity and the spectral type ( Puis et al. 1996 ; Kudritzki 
et al. 1999). Thus observation and analysis of stars of different 
spectral type and metallicity at known distances are crucial to 
calibrate the method before applying it. 

For this reason we have started in our group a number 
of programs aimed at calibrating the WLR for different stars 
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al. 2001). 

In the present paper, we analyze STIS HST UV spectra of 
M33 B-supergiants. Our analysis is primarily aimed at obtain- 
ing the wind terminal velocities, needed to calculate the WLR, 
that will be determined after the analysis of the optical spec- 
trum. 

In addition our study cons titutes, together wi th a similar 
study of B-supergiants in M3 1 ( Bresolin et al. 2002 ), the largest 
homogeneous sample of stars analyzed in the UV beyond the 
Milky Way and the Magellanic Clouds (see Prinja & Crowther 
1998, for a recent list), and an important increase in the to- 
tal number of analyzed B-supergiants, in particular those with 
large intrinsic brightness (bright magnitudes are of course stiU 
a primary selection criteria for our extragalactic targets). The 
group of early B-supergiants is particularly important for the 
spectral type dependence of the WLR. Between 09 and B3, 
two jumps are observed in the WLR, one between 09 and 
BO and another between Bl and B1.5, while the WLR for B- 
supergiants between B1.5 and B3 is not well understood in 
terms of the radiatively driven wind theory, as wind momenta 
l ower than expected a re derived for the analyzed Galactic stars 
dKudritzki et al. 1999|). 



In Sect. ^ we show the observations of the spectra, that are 
described in Sect. Sect. ^ presents the analysis of the U V pro- 
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files, while in Sect. p| we discuss the individual aspects found 
for each star. In SectT^ we present the discussion of the results 
and the conclusions. 



2. The observations 

Observations were made with the HST STIS, using the GL140 
grating which provides a resolution from 310 to 210 km s^^ 

in the wavelength range from 1 150 to 1700 A. We used a 0. 2 
wide slit, which is recommended for optimizing the spectral 
purity. Table |l] gives our list of objects and other details of the 
observations. 

There are two spectra for each star We merged these two 
spectra after checking their relative displacement by cross- 
correlating them. 

Merged stellar spectra are corrected for the relative velocity 
between stars and observer This is done by measuring the dis- 
placement from rest wavelength of metal lines that are purely 
photospheric. There are a few photospheric lines that could in 



principle be used in the observed spectral range (Prinja 1990) 



They are difficult to identify in all spectra. Thus, after inspec- 
tion of the spectra, we decided to use the Si III doublet at AA 
1500.24, 1501.19 A, because it is the only one that can be 
clearly identified in all spectra. The spectral resolution does 
not allow us to resolve both components of the Si III doublet. 
We use as rest wavelength that of a composed line centered at 
A 1500.72 A. Because of the small signal to noise ratio of the 
spectra, we decided to correct one of the spectra of higher qual- 
ity and then refer the others to this one. We selected the spec- 
trum of M33-0900 for that purpose. Once the M33-0900 spec- 
trum had been shifted so that the Si III line was at rest wave- 
lentgh, we checked this correction with the position of other 
strong interstellar lines (see f.e. Prinja 199C ), as these are prob- 
ably produced by the interstellar medium close to M33-0900. 
We can see in Fig. |l]that the shifted interstellar lines are at their 
expected rest wavelengths. 

To correct the other spectra we cross-correlated them with 
the corrected spectrum of M33-0900. To avoid possible biases, 
stron g P-Cygni profiles hav e been masked]^ before correlations 
(see iHowarth et al. 1997|). For M33-1137 and M33-B-526, 



the stars with lowest SN spectrum, the cross-correlation func- 
tion does not present a well defined maximum. Checking the 
strong interstellar lines present in the spectra (Si II A A 1260.40, 
1526.70) we found good concordance between all spectra, ex- 
cluding M33-1137 and M33-B-526. We decided to shift both 
spectra until the position of their interstellar lines agree with 
those of other spectra. Special care was taken in the case of 
M33-B-526, as the unshifted spectrum shows several lines in 
the vicinity of the rest wavelength of the intestellar lines. 

After having shifted the spectra we rectified the continuum 
by tracing a polynomial through a number of selected contin- 
uum points chosen i teratively, in the same way as described in 
Herrero et al. (200l|). 



3. Description of thie spectra 

After all the above corrections have been performed, we finally 
obtained the spectra displayed in Fig. [l|. 

One of the difficulties of stellar analyses in M33 is the large 
possible variation in metallicity, that complicates the interpre- 
tation of the spectrum. Thus in Fig. ^ we show a comparison 
of the N V, Si IV and C IV fines of the M33 stars with Galactic 
stars of similar spectral classification. Their identifications and 
spectral types are given in Table ||. 

The Galactic UV spectra have been taken from the INES 
lUE archive, and have been degraded from the original high 
resolution to the lower resolution of our M33 HST spectra. We 
emphasize that we have taken the rectification points of the 
M33 stars as a guide to rectify the lUE spectra in the same 
way. Therefore, the extra absorption seen redwards of the Si IV 
and C IV profiles in some Galactic stars are probably an effect 
of metallicity. Furthermore these extra absorptions are consis- 
tent with the profiles indicating a lower metallicity in the M33 
stars. The N V comparison is strongly affected by the different 
La absorptions. 

M33-0900 has Si IV and C IV profiles that are very similar 
to those of HD 154090. Although N V is stronger in M33-0900, 
the star looks normal for its spectral type. 

This is not the case for M33-110-A, which in comparison 
to HD 148688 displays very weak Si and C lines. Moreover 
considering the luminosity effect we would expect for Si IV, the 
difference is even larger. This is consistent with the low Si and 
O abundances derived by Monteverde et al. (200C). Thus we 



can expect a low metallicity for 1 10- A, also consistent with the 
higher continuum level redwards from Si IV and C IV. The N V 
profile is again stronger than in the comparison star. The slopes 
of the blue wings of both the Si and C lines are smaller than in 
HD 148688. This will be interpreted during the fit procedure as 
a larger dispersion velocity. Finally we shall mention the broad 
absorption in the C IV blue wing that becomes more apparent 
due to the low absorption in this line. This additional absorption 
is present in both individual M33-110-A spectra and can be 
clearly indentified in Fig. |^ where we show a comparison of 
the C IV line in 1 10-A and in B-38. We have used this figure as 
a guide in fitting the C IV line in 1 10-A. 

M33-B-38, for which [Monteverde et al. (2000|) again de- 



rive low O and Si abundances, has also been compared with 
HD 148688. Similar comments as for 1 10-A are applicable, al- 
though effects are less extreme and the N V profile in this case 
is marginal (as in the Galactic star). 

This has previously been observed and analized by Bianchi 
et al. (1996) and was later reanalyzed by Prinja & Crowther 
(1998). Both used the same HST GHRS spectrum, which looks 
like ours, including the lower slope in the blue wings as com- 
pared to Galactic stars (see Fig. 4 from Bianchi et al. 1996). 
Bianchi et al. find a narrow absorption feature (nar- 

row absorption component, NAC) at a blueshifted velocity of 
^1250 km s^^, which is close to the edge velocity determined 
by Prinja & Crowther (1998), 1225 km s^^. Interestingly, we 



Spectral regions of N V, C IV, Si IV UV resonance lines have been 
set to zero 



also find evidence of a NAC at the same velocity, which may 
indicate that this is a more permanent feature. 
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Table 1. M33 stars observed with HST/STIS. The first column gives the star identification in the catalog of Ivanov et al. 
(1993) (IFM93), while the second one gives the identification by Humphreys & Sandage (1980) (HS80). Coordinates are taken 
from [vanov et al. (1993) and precessed for equinox 2000., except for B-38 and B-133, for which coordinates from Massey et al. 
(1996) are taken. Magnitudes are from the same sources, while spectral types are from Monteverde et al. (1996) and Monteverde 
et al. (2000) , except for 0900 and 1 137, for which new spectral types are given (that will be justified in a forthcoming paper) 



Went 


Went 


a(2000) 


(5(2000) 


V 


Spectral 


Obs. 


Exp. 


S/N 


IFM93 


HS80 






mag. 


Type 


Date 


time(s) 




0900 




01 33 44.7 


30 36 18 


17.3 


BO-Bl I 


Jul-8-00 


5193 


24 


0785 


110- A 


01 33 41.0 


30 22 37 


16.1 


Bl la-^ 


Aug-7-99 


5201 


29 




B-38 


01 33 00.8 


30 35 05 


16.7 


Bl la 


Sep-5-00 


5205 


30 


0515 


B-133 


01 33 29.0 


30 47 44 


17.6 


B1.5 la 


Aug-5-99 


5205 


30 


1733 


B-526 


01 34 15.9 


30 33 46 


17.6 


B2.5I 


Dec- 11-00 


5201 


15 


1137 




01 33 53.1 


30 35 28 


16.7 


B3I 


Jul-9-00 


5205 


15 




-2000 2000 

VELOCITY IK KM/S 



Fig. 3. Comparison of the C IV line in M33-1 10-A (solid line) 
and M33-B-38 (dashed line) 



M33-B-133 shows the most puzzling spectra. The red com- 
ponent of the Si IV doublet is deeper than the blue one, which 
cannot be explained by the theory. The effect is due neither to 
noise nor to variability, as it is present in the same form in both 
individual spectra of M33-B-133. We see numerous features 
in the C IV profile that could be interpreted as NACs. From 
them, only the absorptions at 250, 585 and 1990 km s~^ show 
both components. None of these NACs, however, seems strong 
enough to explain the red component of Si IV. 

This is not the only remarkable feature in the spectrum of 
M33-B-133. This object has been classified by Monteverde et 
al. (2000) as B1.5 la, and is compared with HD 152236 (B1.5 
la-n) in Fig. ^. We see that the profiles are quite different. 
In fact, the spectrum of B-133 resembles more that of n Ori 
(HD 38771), a BO. 5 la star, although a comparison of the opti- 
cal spectra excludes such an early spectral type for M33-B-133. 

Another possibility is a blend with the spectrum of a hot- 
ter star. This is supported by the resemblance of the spectra of 
B-133 and n Ori, and by the data in Table |l| Note that there 
B-133 has a larger S/N than would be expected from a compar- 
ison of its visual magnitude and spectral type with those of the 



other stars. Thus, although a particular set of circumstances is 
required because the B 1 .5 supergiant dominates the blue spec- 
trum, it is not impossible that a hotter companion of type 09.5- 
B0.5 is dominating the UV spectrum. If this is the case the 
apparent strong red component of the Si IV doublet would be 
the main contribution of the B 1 .5 star in the UV spectrum. 

However, we could not reproduce the observed spectrum 
by combining lUE spectra of different spectral types without 
extra ad hoc assumptions (as arbitrary relative displacements), 
nor can we detect any clear indication of binarity in the HST 
WFPC2 images of B-133, or other features indicative of a com- 
panion in the blue and UV spectrum, and thus the suspected 
binary nature of B-133 remains unconfirmed. 

M33-B-526 displays lines very similar to those of HD 
198487, in spite of its suspected binary nature from its appear- 
ance in the Keck I screen (McCarthy, private communication). 

The C IV profile displays several absorption features red- 
wards of the rest wavelength that could be interpreted as con- 
tributions from a companion star. However, their positions do 
not correspond to those of the doublet separation Thus we ex- 
pect this to be the spectrum of a single star. 

M33-1137 has SilV profiles that are clearly weaker than 
those of the selected Galactic counterpart, HD 51309. The C IV 
profile in M33-1 137 is clearly broader than in HD 51309, with 
a shallow slope in the blue wing. It displays additional absorp- 
tions that agree with the positions of the C IV components at 
—75 km s^^, thus clearly pointing to wind inhomogeneities. 
Weaker absorptions are compatible with components at —550 
km s^^. Other absorptions cannot be clearly assigned to the 
C IV doublet, specially considering the low SNR. The N V fea- 
ture is more evident in M33-1 137 than in the Galactic star. The 
observed spectrum points again to a lower metallicity in the 
M33 star than in the Galactic comparison star. 

Therefore the UV spectral morphology of the observed 
M33 stars separates them in three groups: a group (particularly 
evident in Si IV) with Unes weaker than their Galactic compar- 
ison stars (M33-110-A, M33-B-38 and M33-1137), a second 
group with lines similar to their Galactic counterparts (M33- 
0900 and M33-B-526) and one star (M33-B-133) showing pe- 
culiar profiles for its spectral type. In addition, the M33 stars 
seem to have lower slopes in the blue wings, pointing to larger 



4 



Urbaneja et al.: STIS UV spectroscopy of M33 early B supergiants 



0) 




1200 



1300 



1400 

Lambda (A) 



1500 



1600 



1700 



Fig. 1. The observed rectified UV spectra. From top to bottom: M33-0900, M33-1 10-A, M33-B-38, M33-B-133, M33-B-526 and 
M33-1 137. Relative fluxes have been arbitrarily displaced in ordinates for the sake of clarity. The main IS lines are marked at the 
top, and below we have indicated the rest wavelengths of the most important stellar lines. 



dispersion velocities. This is particularly clear in the first group 
of stars. 



4. Analysis of the resonance lines 



We use the method described by Haser (Haser 1995 , see also 



Lamers et al. 1999) to analyse the UV P Cygni lines in order to 
derive wind terminal velocities. We have also derived H I col- 
umn densities towards the M33 stars by fitt ing the IS La line, 
in the same way as i n ^r[errero et al. (2001 ) (see also Jenkins 
1970; Bohlin 1975 ). The values quoted in Table ^ for these 
column densities are similar to those obtained for stars in our 



Galaxy (see Fig.2 of |Shull & Van Steenberg 1985D . 

After continuum rectification, we extract small regions of 
each spectrum centered in our strategics doublets: Nv AA 
1238.819, 1242.798; Siiv AA 1393.73, 1402.73 and Civ AA 
1548.191, 1550.761. In case of discrepancies or difficulties, we 



give more weight to the Si IV and C IV lines than to the N V 
one. 

To get the wind terminal velocities we have to correct for 
the underlying photospheric components, which we do in an 
approximate way, by using lUE spectra of hot stars with weak 
winds (and projected rotational velocities as low as possible) as 
templates. These templates have been convolved with appro- 
priate rotational profiles, to account for the individual stellar 
rotational speeds. We excluded N V from this procedure, as the 
corresponding spectra are noisy and do not show any features. 

We selected a sample of Milky Way dwarfs with spectral 
types from 09V to B3V taken from the INES database. This 
covers the spectral type range of our M33 stars. Selection of 
photospheric template and continuum rectification has a big 
impact when fitting emission peaks, but has little effect on the 
determined terminal velocities. The lUE spectra have been cor- 
rected for the relative velocity between stars and observer, in 
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Fig. 2. Comparison of the Nv (left), Si IV (middle) and CiV (right) lines in the M33 (full Une) and Galactic (dashed line) 
stars. From top to bottom: M33-0900 and HD 154090; M33-110-A and HD 148688; M33-B-38 and HD 148688; M33-B-133 
and HD 38711; M33-B-526 and HD 198487; M33-1137 and HD 51309. Spectral clssifications are given in Table |. Relative 
fluxes have been arbitrarily displaced in ordinates for the sake of clarity. Galactic spectra have been degraded to the M33 HST 
resolution. 



the same way as for M33 stars. The stars used for the photo- 
spheric templates are specified in Table ^. Note that the profiles 
of interest are very weak in these stars, but even so the different 
metallicity may be a small source of error. 

The velocity stratification is parameterized as a usual /3- 
law, and we account for the wind turbulent velocity in the same 
manner as described in Haser (1995) and Herrero et al. (2001 ). 



5. Individual comments 



We briefly comment here on individual aspects of the analysis 
that could be of interest. 



The indetermination in the exponent of the velocity field, f3, 5.1 . M33-0900 
produces an additional uncertainty in the terminal velocities. 

Table ^ gives the fit results with respect to the wind and 
turbulent velocities. The final fit to each line is shown in Figs. Q 
and I 



We used HD 22591 (B0.5 V) for the photospheric profiles. The 
N V profile shows a large emission peak that we cannot repro- 
duce completely (see Sect. ||). The Vta/voo ratio is 0.18, larger 
than in Galactic stars. 
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1,2 




Lambda {A; Lambda (A) Lambaa (A) 

Fig. 4. Final fits for the three earliest stars of our sample. Lines are plotted from left to right (N V, Si IV, C IV) and stars from top 
to bottom in the same order as they are listed in the tables. 



5.2. M33-110-A 

We used HD 39777 (B1.5 V) for the underlying photospheric 
profiles. We find the same problems as for M33-0900: a large 
red emission peak in N V, and a large ratio vta /voo, which with 
a value of 0.25 is much larger than in Galactic B-supergiants. 
Furthermore, the additional absorption contaminating the blue 
side of C IV makes the fit a bit more difficult, affecting mainly 
the determination of Vta. This does not seriuosly affect the un- 
certainty in the adopted values for (3, v^o and Vta- 



5.3. M33-B-38 

The terminal velocity of this star has also been determined by 
Bianchi et al. (1996| ) and ^rinja & Crowther (1998"| ). The first 
authors use an analysis method similar to the one employed 



here, and thus it is not surprising that our value agrees with 



their. Prinja & Crowther (1998) determine the terminal veloci- 



ties of their stellar sample from the violet edges of the profiles 
and the NACs. The velocity they obtain for this star is much 
larger (1225 km s~^) than ours (730 km s~^) or the one by 
Bianchi et al. (about 700 km s~^). Our velocity would support 
the interpretation that what we have seen at 1250 km s^^ is 
actually not a NAC. 



Bianchi et al. (1996| ) do not give the turbulent velocity of 
their stars. We derive a very large value, 250 km s^^, nearly 
35% of Voo- 

We adopted a large /9-value, both to improve the consis- 
tency between the Si IV and C IV fits, and to improve the fit of 
the blue side of the red emission peak. 
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Fig. 5. As Fig. ^ however for the three last stars in Table |l] 



5.4. M33-B-133 

We have used HD39777 for phostospheric profiles, as in the 
case ofM33-110-Aand M33-B-38. 

In spite of the low O and Si abundances derived by 



Monteverde et al. (200C ) for this star, its UV spectrum shows 
comparatively strong C IV and Si IV P-Cygni profiles, indicat- 
ing a strong and fast wind (see Sect. |]). The terminal velocity 
reaches 2050 km s~^. The turbulent velocity is 150 km s^^, a 
modest 7%. This is the only star for which we obtain a terminal 
velocity clearly above the galactic average for its spectral type 



(taken from Kudritzki & Puis 200C ), thus challenging the low 



abundances or the spectral classification (or both!). The spec- 
tral type, however, should change from B 1 .5 to 09 in order to 
have a terminal velocity close to the spectral type average. This 
is completely ruled out from inspection of the optical spectrum. 



The fits displayed in Fig. ^ also show a behaviour differ- 
ent from those of the other stars. N V is well fitted. The fit of 
the C IV profile is good in the bluest part of the wind absorp- 
tion profile but is bad in the rest of the profile. Fortunately, the 
first one is the important part for determining the terminal ve- 
locity. The unsatisfactory fit at low wind velocities is cosmet- 
ically very dependent on the underlying photospheric profile, 
and thus does not mean very much by itself, but it would be 
consistent with the presence of additional C IV absorption at 
low wind velocities. We keep the low /3 value as there is no 
need to increase (3 to improve the consistency between C IV 
and Si IV. However, the anomaly in Si IV indicated in the pre- 
ceding section, i.e., a stronger red component, can also not be 
fitted and indicates an extra absorption at low wind velocities. 
This uncertainty does not affect the determination of the ter- 



8 



Urbaneja et al.: STIS UV spectroscopy of M33 early B supergiants 



minal and turbulent wind velocities, that have an accuracy of 
±100kms-i. 

This does not solve the problem referred to above. The star 
has a large v^c for its spectral classification and a Si IV red com- 
ponent stronger than the blue one. The first might be attributed 
to a case of bi-stability ( [Pauldrach & Puis 199C ), similar to that 
of P-Cygni in our Galaxy, where the ratio of terminal velocity 
to escape velocity is larger. This however does not explain the 
red Si IV component. One possible solution to this puzzle is that 
we are looking at a composite spectrum, B-133 being the star 
contributing to the red component of Si IV (see Sect. ||). If this 
is the case, we estimate its terminal wind velocity to be of the 
order of 450 km s~^, but would expect to see extra absorption 
in Si IV at low wind velocities. 

5.5. M33-B-526 



It has been classified as B2.5I ( [VIonteverde et al. 1996 ), from 
its optical spectrum. We used a B2.5 V galactic star for the pho- 
tospheric profiles (HD 44402, Z CMa). 

The derived terminal velocity is the lowest in our sample 
and the UV spectrum shows only weak signs of mass loss. 
However, the turbulence velocity is relatively high, with a ra- 
tio vta/voo in excess of 0.3. This large value is the result of a 
compromise between the fit of C IV and that of Si IV. The first 
would allow a lower turbulence, but then we cannot fit the sec- 
ond one. Again, the value of (3 has to be large to favour consis- 
tency between both doublets and a better fit to C IV. However, 
the terminal velocity is not seriously affected. The uncertainty 
of both velocities is different now, being that of vta ±100 km 
s~^ and that of Voo ±50 km s~^. For the sum of both we have 
adopted ±75 km s^^. 

5.6. M33-1137 

This is the coolest star of our sample, with spectral type B3Ia. 
We take HD32630 (B3V) for the photospheric profiles. 

The Si IV doublet is mainly photospheric or in any case the 
wind has a small contribution. The fit to N V is poor. 

The C IV is the main profile for deriving the parameters of 
the velocity field. The low SNR of the spectra and the addi- 
tional absorptions described in the preceding section make the 
fit difficult. We find the best fit at 750 km s"\ 250 km s^^ and 
2.0 for Voo, Vta and /3. Again, Voc is larger than the Galactic av- 
erage from Kudritzki & Puis (2000), but now much more mod- 
erately than for B-133. The uncertainties are slightly larger than 
those of the other fits, as akeady expected from the described 
difficulties. We adopt ±100 km s~^ for v^o, Vta and its sum. 



6. Discussion and conclusions 

In Fig. ^ we have plotted the derived Voc against the stellar 
spectral types, together with parabolic fits to the average val- 



ues quoted by Kudritzki & Puis (2000) for Galactic OB su 



pergiants. The fits have been obtained by joining 01 and Bla 
supergiants on the one hand and Oil and Bib supergiants on 
the other. We have also plotted lines that indicate a 30% vari- 




Fig. 6. The Voo obtained for the M33 stars. Solid lines repre- 
sent the fits to 01 + Bla (lo wer curve) and to Oil + B ib super- 
giants (upper curve) of the Kudritzki & Puis (2000| ) Galactic 
data. Dashed lines are ±30% variations of these curves. The 
solid dot indicates the approximate position of B-133 if derive 
it from fitting the red Si IV component (see text for details) 



ation from th e plotted average relati ons, a usual scatter range 
(see Fig.4 in [Kudritzki & Puis 2000| ). 

All M33 B-supergiants give us values that can be consid- 
ered normal, except B-133. We see no difference between the 
stars with different suspect ed metallic ities. This is in agree- 
ment with [Puis et al. (200C|) (see also |Vink et al. 200C|) who 



have argued that the terminal velocity (which depends on the 
slope of the line-strength distribution function, a) is primarily 
controlled by the ratio of light ions vs. iron group elements, 
because of the different line statistics. As long as this ratio is 
roughly similar, the theoretical expected change of terminal ve- 
locity (due to the "indirect" a effect, see Puis et al., 2000, setc 
5.2) is much smaller than if this ratio would be changed. In 
particular if a dense wind is present (as is the case for our su- 
pergiants) and this ratio remains unchanged (which we have 
to assume for the moment), the effect is expected to be very 



small (see Puis et al. 2000, Fig. 27), since the effective a then 



remains roughly constant. Only for thin winds and/or a sig- 
nificantly lower (general) metallicity the effect should become 
observable. 

We could not find a satisfactory explanation for B-133. 
While its terminal velocity is well above the average for its 
spectral type, it could still be accepted (even with the low 



metallicity derived by Monteverde et al. 200C) assuming for 



example a bi-stability phenomenon, were it not for the anoma- 
lous red component of the Si IV doublet. Assuming on the other 
hand that B-133 is actually producing the red Si IV component, 
its velocity would be much closer to the average of its spectral 
type. However, we could not reproduce the observed profiles 
without ad hoc hypothesis, nor find conclusive evidence of a 
binary nature. 

Although three of our stars have been analyzed by 



Monteverde et al. (2000) we will postpone a discussion of the 
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Table 2. Results obtained for the observed stars. Voo and vta (turbulence velocity in outer wind) and the uncertainty of their sum, 
Av, are given in km s~^ 



IQCIIL 








log 






Z_i u 


P 




Type 


Star 


Type 


N(Hl) 










M33-0900 


BO-Bl I 


HD 154090 


B0.7 la 


20.85 


950 


170 


50 


1.0 


M33-110-A 


Bl Ia+ 


HD 148688 


Bl la 


21.15 


800 


200 


50 


1.0 


M33-B-38 


Bl la 


HD 148688 


Bl la 


20.70 


730 


250 


50 


2.0 


M33-B-133 


B1.5 la 


HD 152236 


B1.5 Ia+ 


20.97 


2050 


150 


100 


1.0 


M33-B-526 


B2.5I 


HD 198487 


B2.5Ia 


21.10 


380 


120 


75 


2.0 


M33-B-1137 


B3 la 


HD51309 


B3Ib 


20.85 


750 


250 


100 


2.0 



terminal velocities in terms of the stellar parameters, as new 
analyses of their optical spectra, now including mass-loss ef- 
fects, are currently under way in our group together with a set 
of newly observed stars. 

The wind turbulent velocities derived here are larger than 
the typical 10% (|Kudritzki & Puis 2000| ) or the 14% found 
by ^errero et al. (2001 ), reaching nearly 35% in the most 
extreme case. Including B-133 we obtain 0.25 for the mean 
Vta/ Voo ratio. This result is of the same order of that found by 



Bresolin et al. (2002) in their analysis of M31 B supergiants. 



We tentatively attribute it for the moment to the extreme char- 
acter of our objects, selected among the brightest supergiants 
in M33, and as a consequence are among the most luminous 
objects. 

Finally, we have detected evidence of numerous NACs, 
confirming the wide presence of wind inhomogenities. 
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